Negativity and quantum discord in Davies environments 
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We investigate the time evolution of negativity and quantum discord for a pair of non-interacting 
qubits with one being weakly coupled to a decohering Davies-type Markovian environment. At 
initial time of preparation, the qubits are prepared in one of the maximally entangled pure Bell 
states. In the limiting case of pure decoherence (i.e. pure dephasing), both, the quantum discord 
and negativity decay to zero in the long time limit. In presence of a manifest dissipative dynamics, 
the entanglement negativity undergoes a sudden death at finite time while the quantum discord 
relaxes continuously to zero with increasing time. We find that in dephasing environments the 
decay of the negativity is more propitious with increasing time; in contrast, the evolving decay of 
the quantum discord proceeds weaker for dissipative environments. Particularly, the slowest decay 
of the quantum discord emerges when the energy relaxation time matches the dephasing time. 
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I. INTRODUCTION 

Establishing correlations is a sine qua non condition in 
effective communication. If two parties are quantum cor- 
related one can attempt to conduct a quantum commu- 
nication between them. Quantum communication proto- 
cols [l| make use of certain properties of states of quan- 
tum composite systems. The best known and most popu- 
lar resource for quantum communication is quantum en- 
tanglement 0, a widely studied property of composite 
systems. Quantum entanglement 0], interesting per se, 
has been recognized as a powerful tool for quantum infor- 
mation processing [l|. For example, let us mention the 
use of quantum communication with dense coding and 
teleportation as most celebrated examples 0, [|[ . 

If one presupposes that the history of quantum en- 
tanglement started out with the work by Schrddinger 
Q, it took almost a century to discover that there are 
quantum correlations which are essentially non-classical 
and these can exist even in the absence of entanglement 
[HHl- When compared to entanglement, quantum corre- 
lations are, to some extent, 'more mysterious', since their 
mathematical setting has not been uniquely established. 
The well known mathematical language for entanglement 
incorporating tensor products of state spaces, Schmidt 
decomposition, etc., not only allows one to pose many 
of fundamental problems in purely algebraic context but 
also motivates the search for a variety of entanglement 
quantifiers. The physical context of entanglement en- 
coded e.g. in Bell-type inequalities Q can be (sometimes 
incorrectly) considered as a consequence of mathematical 
theorems. Recent studies of quantum correlations have 
shown that a broad class of composite systems carry cor- 
relations which can be described in the context of infor- 
mation extracted from composite parties via a suitably 
formalized measurement procedure, see in Ref. 1 for re- 
cent review. 

Among several quantifiers of quantum correlations the 
quantum discord attained recent popularity^. Despite its 
fundamental meaning, exemplarily studied via approxi- 



mating a reduced quantum dynamics within trace pre- 
serving completely positive maps [6], quantum correla- 
tions quantified by quantum discord can open new av- 
enues for quantum computations and quantum com- 
munication schemes [8l-[l0j|. 

A salient obstruction in the implementation of both 
quantum computation and communication protocols is 
quantum decoherence [l| , induced by the influence of the 
environment on the relevant 'quantum hardware'. Stud- 
ies of entanglement dynamics in the presence of noise 
enjoy a long history [ll|, Qjj]; in contrast, similar studies 
on quantum discord dynamics have been carried out only 
recently [13l4la |. 

According to Ref. Qjl, the quantum discord and the 
entanglement present fundamentally different resources. 
Here, we investigate if this difference manifests itself 
in their robustness to environment-induced decoherence. 
In other words, we study whether it is possible to as- 
sign an environment to two general classes, one of them 
more suitable for quantum information processing using 
entanglement -based protocols, and the second preferring 
quantum correlations as being quantified by quantum 
discord. With this study we obtain an affirmative an- 
swer to this question. 

In doing so we consider a simple set-up: We investi- 
gate the time evolution of quantum correlations between 
a pair of qubits with only one of them being coupled to a 
decohering environment. We emphasize here that the ex- 
act reduced quantum dynamics for the two qubits is typ- 
ically neither comp letely positive 0Qil, nor is it gener- 
ally even linear [20| , not to speak of undergoing a memo- 
ryless quantum Markovian dynamics [21] . Here, however, 
we restrict ourselves to a quantum Markovian dynamics, 
having in mind quantum optics with weak coupling at 
not too low temperatures. Then the time evolution can 
satisfactorily be approximated with a trace preserving 
completely positive Markovian map of Davies-type [22J . 

The structure of the work is as follows: In Sec. 2, we 
describe our model set-up and its dynamics in terms of 
a Davies map. Next, in Sec. 3, we consider the time 
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evolution of the correlation quantifier negativity while 
with Sec. 4, we analyze the quantum discord. With Sec. 
5 we present our summary and conclusions. 

II. OPEN QUANTUM DYNAMICS: QUANTUM 
MARKOVIAN DAVIES MAP 

With this work, somewhat stimulated by prior work 
done in Ref. fl3l . we consider a pair of non-interacting 
qubits A and B that initially are prepared in one of the 
maximally entangled pure Bell states; i.e., 

Pi = (vo®Vi)Po(vo®Vi),i = 0,1,2,3 (1) 

where 

Po = ~(|01) + |10))(<01| + <10|). (2) 

The non-decohering part of the evolution of the qubits A 
and B is determined by the two two-level Hamiltonians; 
i.e., 




For the following, only qubit B is coupled at times t > 
to an environment E at temperature T . Then, in the 
presence of this environment coupling, the dynamics of 
the two qubits, described by the reduced density ma- 
trix pab (t) , no longer proceeds unitary. Because the two 
qubits A and B are initially correlated, both the quan- 
tum entanglement and the quantum discord will evolve in 
course of evolving time t > 0, as these become influenced 
by a finite system B-environment coupling. 

In presence of a very weak system-reservoir interac- 
tion and not extremely low temperatures the deduced 
dynamics can satisfactorily be described [2l| by a Marko- 
vian dynamics of the Davies type [22| . The main advan- 
tage of this approach is that it recovers the well estab- 
lished steady state properties at assumed weak coupling, 
namely stationarity and asymptotic stability in terms of 
the Gibbs state, reading 

Pfl(*)-> ({J ilp) for oo. (4) 

Here, the thermal weight is p = exp(— f3u>B)/Z, Z — 
exp(— (3lub) + exp(/3ws) and /3 = 1/kT denotes the in- 
verse temperature. Using Davies theory one explicitly 
can construct the generator of a completely positive semi- 
group describing the reduced dynamics (with respect to 
the environment) of the open quantum system in terms 
of microscopic parameters in the Caldeira-Leggett-type 
Hamiltonian of the full system. Such modeling has re- 
cently been applied to studies of entanglement dynam- 
ics [23|], properties of a dissipative geometric phases of 
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FIG. 1: (Color online) The sketch of the system set-up studied 
with this work. The qubits A and B with level separation uja 
and lob, respectively, do physically not interact; these are, 
however, initially prepared in one of the maximally entangled 
Bell states, as indicated by the broken (green) line. Here, 
only the qubit B is weakly coupled to an environment. The 
dynamics of the qubits is given by Eq. ((5]) in the text. 



qubits ^24] and for some thermodynamic properties of 
nano-scale systems [25j . to name but a few. 

Here, instead of using the most general scenario of 
Davies semigroups, we restrict ourselves to a completely 
positive Davies map [Hj]. The qubit-qubit reduced dy- 
namics is given by the following Davies map $f [2q . 
acting on the Hilbert space B and completed with the 
tensorized unitary time evolution for the uncoupled qubit 
A, the latter acting solely on the Hilbert space of qubit 
A. Put differently, we have for the reduced dynamics a 
(super-operator)-time evolution, reading 

PAB(t)=tpt®*t)pi, » = 0,1,2,3. (5) 

This is because the A-system Hamiltonian Ha commutes 
with all other remaining Hamiltonians of the system B 
and environment E. Stated explicitely, for any linear 
combination of matrices the form <8> |*b)0'sI> we 

have (U?®$?)\i A )(j A \ ® \i' B ){j' B \ = (|*a>0'aI) ® 
3>f (Ks)Osl)- The part of the super-operator U* is the 
unitary evolution operator generated by the Hamiltonian 
H A ; i.e., U t (\i A )(j A \) = e- iH ^\i A ){j A \e iH ^. Let us de- 
scribe how the non-unitary Davies map <FP acts on the 
one-qubit density matrix p B - Following |26| let us con- 
struct a super-operator &f corresponding to <&f acting 
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on vectorized matrices 



Pb = 



Poo Poi 

PlO Pll 



\pb)) = 




$f p B -► #f \\ PB )). 



(6) 



(7) 



An explicit form of the Davies super-operator then reads 
261: 



l-u{t) u(t)p/(l-p) 

v(t) 

v*(t) 

a 1 - u(t)p/(l - p) 



(8) 



with the corresponding relaxation functions reading 

u{t) = (l-p)[l-exp(-Ft)], 
= cxp(—iu!Bt — Gt). 



(9) 
(10) 



The parameters i* 1 = 1/ri and G = 1/t2 are related to 
the energy relaxation time T\ and the dephasing time Ti , 
respectively. Given the fact that these relaxation times 
are subjected to obey the physical inequalities 27] 



G > F/2 > 0, i.e. 2ri > r 2 , 



(11) 



it then follows that this Davies map indeed is a trace- 
preserving, completely positive map. 

Let us consider the limiting case of an infinite energy 
relaxation rate, yielding F — 0. This corresponds to pure 
dephasing without relaxation in energy taking place (no 
dissipation). Such dephasing scenarios, including also 
non-Markovian dephasin g m odels, have been applied to 
studies of entanglement |28j ] and quantum discord [29| 
dynamics. The opposite case, i.e. G = and F ^ 
can physically not be realized: If there is dissipation of 
energy then necessarily finite dephasing accompanies this 
relaxation process. For p = 1/2 in Eq.(|4]), i.e. in the 
case of an infinite large temperature, the single qubit 
B state becomes maximally mixed: Tr A(pAB(t)) = |Ib- 
Likewise, we have for qubit A that = TrB(pAB(t)) — 
kIa. Hence with p — 1/2 the quantum discord, being 
evaluated below, is symmetric with respect to A and B 
labelling 

Because the explicit results below are robust with re- 
spect to any chosen value for the Boltzmann weight p, 
we shall depict in our Figures the case with p = 1/2 only, 
yielding a symmetric quantum discord. 

Given that at initial time the reduced two-qubit den- 
sity matrix pab{0) = Po, where po is the Bell state given 
by Eq. the density matrix pab (t) evolving under 

the Davies map |5J assumes the general time-dependent 
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FIG. 2: (Color online) Upper panel: The entanglement quan- 
tifier for negativity N(t) [see Eq. (|14[l ] at temperature T = oo, 
i.e. p — 1/2, as a function of the dimensionless time Gt for 
several values of the ratio a = F/G = ts/ti, where ti and 
n are the dephasing and the energy relation times, respec- 
tively. The inset depicts the non-exponential decay of nega- 
tivity N(t) for a dissipative environment, implying that F > 
(yielding a > 0). The negativity N(t) vs. the dimensionless 
time Gt is depicted on the semi-logarithimic plot with the 
vertical axis scaled logarithmically. Lower panel: Time evo- 
lution of the quantum discord D(t), [see Eqs. (1221) and (1231) ]. 
The inset shows the difference between discord in the dissipa- 
tive case(i.e. F > 0) and for pure dephasing (i.e. F = 0) for 
several values of the ratio a. The slowest decay of quantum 
discord D(t) is observed when ri = Ti. 



main result: 

pAB(t) 
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1 - e~ Ft 



(12) 



where to — uja — and Pab (t = 0) = po- This density 
matrix assumes the form of a so termed X-state [l7l[30|. 
Note that the X-structure of the reduced density ma- 
trix (|12[) is preserved during time-evolution. This feature 
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originates both from the symmetry of initial preparation 
(P), here assumed to be given by Eq.®, and the charac- 
ter of dynamics given by completely positive Davies map. 
Let us stress, however, that none of the physical results 
reported below depend on the choice of the specifically 
chosen initial Bell state in Eq. ([2]); put differently, the 
results remain robust for any of the four Bell states. 

III. ENTANGLEMENT DYNAMICS 

In this section we investigate entanglement of the two- 
qubit state pab (t) given in ([5]). As a proper measure of 
entanglement we use the quantifier of negativity N(t) 
being defined as 

N ® = ^E(i A >i- A *) . ( 13 ) 

i 

where the Xi are the eigenvalues of the partially trans- 
posed density matrix pAB{t), at fixed time t, of a com- 
posite system Q. 

Let us next evaluate this negativity of the state (|T2"j). 
It explicitly reads: 

8N(t) = 2er Gt + er Ft - 1 + |2e~ Gt + e~ Ft - l| . (14) 

N(t) does not depend on to, i.e., it is not affected by the 
individual single-qubit level spacings in Eq. ([3]). The 
time evolution of the negativity N(i) is depicted with the 
upper panel in Fig. [5]for T — > oo. For pure dephasing; i.e. 
for F = the negativity exhibits a strictly exponential 
decay, reading: 

N(t) = ±e- Gt , (15) 

with the characteristic dephasing rate G = I/V2. For a 
dissipative environment with F > 0, the entanglement 
undergoes a sudden death (ESD) [l2| occurring at finite 
death time t c . Put differently, for time t > t c , the entan- 
glement vanished identically. From Eq. (fl4l) it follows 
that t c is determined by the relation 

2e~ Gt " + e~ Ft " -1 = 0. (16) 

Upon increasing the relaxation rate F the death time 
t c monotonically decreases from t c — 00 for F = to 
the minimal value t™ in = In (y/2+l) /G, which occurs 
when F/G = t^/ti = 2. In the case of finite dissipation, 
F > 0, the negativity JV(i)-decay proceeds faster than 
exponentially, see the inset in the upper panel in Fig. 2, 
where N(t) vs. the dimensionless time Gt is depicted on 
the semi-logarithimic plot with the vertical axis scaled 
logarithmically. 

IV. QUANTUM DISCORD 

Next let us investigate a quantum correlation measure 
as encoded with the quantum discord. We again con- 
sider a composite system consisting of the two qubits A 



and B. The full (classical and quantum) correlations in 
the composite system are encoded with quantum mutual 
information, defined as 

Apab) = S{ PA ) + S(p B ) - S(pab), (17) 

where S(p) = — Tr(plnp) denotes the von-Neumann in- 
formation entropy, pab is the density operator of the 
composite bipartite system AB. The part pa refers to 
the reduced density operator of system part A while, like- 
wise, pb is the reduced density operator for system part 
B. 

The classical part of the total correlations is defined as 
a maximum information about one subsystem A that can 
be obtained by performing a measurement on the other 
subsystem B, as defined by a complete set of projectors 
{Ilf }. Let us recall that for p = 1/2, see in Sec. HI an in- 
dependence of the chosen A, B-labelling is granted. The 
label k distinguishes different outcomes of this measure- 
ment. The quantifier of the classical part of correlations 
is defined by the set of following relations 

C( PAB ) = maxS(p AB |{nf}), (18) 
{nf} 

1 

S( PAB \{ILZ}) = '£ qk S(p k A ), 

P \ = — Tr B [(i A ® nf ) pab (Ia ® nf )] , 

Qk 

qk=Tr AB [(l A ®n%)p AB ], (19) 

The difference between the total amount of correlation 
and the classical part of correlation thus reads 

D{ PAB ) =1{pab)-C{ PA b) ■ (20) 

This relation defines the quantum discord. Being so, it 
provides a measure for manifest quantum correlations [5j . 
Let us remark that, generally, this quantum discord (|20[) 
presents neither a unique nor the most optimal quanti- 
fier for quantum correlation Q. However, for the case 
of bipartite systems one can summarize that the states 
can be divided into two groups [6]; namely entangled 
(quantum correlated) and separable states. In turn, the 
separable states can either be classically correlated or 
quantum correlated (but not entangled). This classifi- 
cation is non-trivial, as e.g. classically correlated states 
always lead to completely positive maps while states with 
quantum correlations may give rise to non-completely 
positive maps Q. For our set-up here here, classically 
correlated states exhibit zero quantum discord D(t) = 0, 
while states with quantum correlations exhibit a non-zero 
discord D(t)>0. 

There occurs a natural computational difficulty in eval- 
uating the quantum discord, stemming from the maxi- 
mization procedure in Eq. (|18[) . Fortunately, this task 
becomes feasible for two qubit systems: it is sufficient to 
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consider projective measurements of the form [5j, l3l| 

B ( cos 2 (0/2) sin(0)exp(^)\ 
^ sin(0) exp(-#) sin 2 (6»/2) )' 

nf = i B -nf, (21) 

where {0, 0} is a standard parameterization of a single 
qubit Bloch sphere. This simplification is helpful in con- 
sidering the quantum discord for models which can be 
effectively described in terms of two qubits [HI, [l|| [l?} ■ 

The calculation of quantum discord requires, in gen- 
eral, a careful optimization with respect to (being suffi- 
cient for us) projective measurements (|21|) . Fortunately, 
here the problem is even more tractable. Due to symme- 
try of the system, finding the optimum in (118[) does not 
involve the non-trivial ^-dependence, i.e. the measur- 
ing process ([21]) can be simplified to a single-parameter 
family of projectors. The quantum discord therefore is 4>- 
independent. In addition it exhibits extrema, depending 
of the value of the dissipation parameter F, For F < G, 
the maximum in Eq. (fT5]) occurs at 9 = in Eq. (j2"T|) . The 
quantum discord D(t) consequently reads: 

4D(t) = (e-^ + l + 2e- Gt )ln(e- Ft + 1 + 2 e - Gt ) 
+ (e~ Ft + 1 - 2e- Gt ) \n{e- Ft + 1 - 2e~ Gt ) 
- 2(1 + e- Ft ) ln(l + e- Ft ), F < G . (22) 

To optimize quantum discord within the required regime 
G < F < 2G (see in (TTTI) ') in turn implies for consistency 
that 9 = tt/2 in Eq.([T8]). The final result for D(t) thus 
reads 

4D(t) = (e- Ft + 1 + 2e- G *) ]n(e~ Ft + 1 + 2e~ G *) 
+ 2(1 - e~ Ft ) ln(l - e~ Ft ) - 2(1 - e~ Gt ) ln(l - e~ Gt ) 
+ ( e - F * + 1 - 2e~ G ') ln(e- F * + 1 - 2 e - G ') 
- 2(l + e~ Gt )ln(l + e- Gt ), G<F<2G. (23) 

It follows from Eqs. and (g3J) that the discord D{t) 
monotonically approaches as t — > oo and stays positive 
for finite time, see the lower panel in Fig. 2. Contrary to 
negativity N(t), which is most robust in the absence of 
energy relaxation F — 0, the slowest decay relaxation of 
quantum discord is obtained for the case when F = G, 
i.e. when the two characteristic relaxation times match, 
i.e., T\ — t%. To clarify this behavior let us notice that for 
pure dephasing F — classical correlations, (|18p remain 
constant without further decay; - put differently, Marko- 
vian dephasing does not affect the dynamics of classical 
correlations. 

The difference between the quantum discord and the 
entanglement becomes best visible in the asymptotic long 
time regime when t 3> T\ or t 3> T<i. Here, one obtains 
the limiting behavior 

D(t) ~ 1 [e- 2Gt + e~ 2Ft ] , for G < F < 2G . (24) 



Therefore, for t — > oo the decay rate of quantum discord 
is determined solely by the dephasing time 1/G = ti and 
is independent of the energy relaxation time n . However, 
the entanglement N(t) vanishes identically in this time- 
regime if finite energy relaxation with F > is at work. 
On the other hand, for a pure dephasing F = one finds 
the asymptotic relation 

D{t) ~ e- 2Gt /2 = 2N 2 (t) . (25) 

Although this very kind of relation between quantum 
discord D(t ) and entanglement N(t) likely may not be 
generic [lj, [13], it nicely illustrates the critical role of 
an environment for the asymptotic behaviors of both 
quantifiers. Notably, the ratio N(t)/D(t) exhibits for 
t — > oo a divergence for a pure dephasing, while it 
vanishes identically in the presence of finite energy 
relaxation. 



V. CONCLUSIONS 

With this work we have investigated the time evolu- 
tion of entanglement negativity N(t) and quantum dis- 
cord D(t) for a pair of qubits, with one qubit (B) weakly 
coupled to a decohereing environment. The decoherence 
dynamics for this sub-system has been approximated by 
a Markovian, completely positive Davies semi-group dy- 
namics. We identified two classes of environments that 
impact differently the quantum correlations as quantified 
by entanglement and quantum discord: (i) In the limiting 
case of the pure decoherence (i.e. strict dephasing) both, 
the quantum discord D(t) and the negativity N{t) decay 
exponentially towards zero in the asymptotic long time 
regime. Moreover, there exists an appealing functional 
relation between these two measures, being detailed with 
([23| . (ii) In the case of a dissipative dynamics, the en- 
tanglement undergoes a sudden death at a finite time 
t c while the quantum discord D(t) smoothly relaxes to- 
wards zero at long times. The slowest decay of the quan- 
tum discord occurs when the energy relaxation time T\ 
matches the dephasing time T2 . Our findings may serve as 
a potential guideline for the implementation of quantum 
information, e.g. communication protocols. We have elu- 
cidated which of the two types of correlation measures, 
namely quantum entanglement N(t) or the one encoded 
by quantum discord D(t), can provide an advantageous 
and/or more suitable quantifier for quantum information 
processing occurring in an open quantum systems under- 
going an ubiquitous decoherence dynamics. 
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